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Abstract 

Here we analyze the relationship between quantum contextuality and decoherence in interference 
experiments with matter particles by means of a simple reduced quantum-trajectory model, which 
attempts to simulate the behavior of the projections of multi-dimensional, system-plus-environment 
Bohmian trajectories onto the subspace of the reduced system. This model allows us to understand 
the crossing of the subsystem trajectories as a combined effect of interference quenching and erasure 
of "which-way" information, which can be of utility to interpret decoherence effects in many- 
dimensional systems where full Bohmian treatments become prohibitive computationally. 
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I. INTRODUCTION 



Quantum physics is characterized by striking properties which puzzle and challenge our 
understanding of the physical world. One of them is contextuality, the unavoidable de- 
pendence of the description of a system on the experimental (or contextual) setup, which 
strongly determines the wave function describing the quantum state of that system [l| . This 
strongly contrasts with the behavior displayed by classical systems, whose dynamics is only 

externa 



governed by the eventua 
matter particles Q, S fl, 3, E, H, Is 



forces. The interference of separate diffracted beams of 
constitutes a nice scenario where this property becomes 
very apparent. Consider, for instance, a two-slit diffraction experiment. If no "which- way" 
information is demanded, the diffracted projectiles will distribute beyond the slits at a time 
t according to the probability density 

p t (r) = |^, t (r) + ^ 2 , t (r)| 2 = p M (r) + p 2>t (r) + p int , t (r). (1) 

Here, ipi,t( r ) (i — 1 ; 2) is the wave function describing the particle after only passing through 
slit i and Pi,t(r) = |^j )t (r)| 2 is the probability corresponding to finding that particle at r. The 
particular fringe-like features observed in the interference pattern arise from the interference 
term p int t (r), which is associated with the contextual or experimental information "both slits 
are open simultaneously when the particle beam reaches them". Otherwise, if we choose our 
context to be one slit open at a time, particles will distribute according to the classical-like 
probability distribution 

p t (r) = |^i, t (r)| 2 + |^ 2 ,(r)| 2 = p M (r) + p 2 , t (r), (2) 

which obviously does not contain the interference term. Thus, for the same physical system, 
two different contexts provide two different answers — note that in classical mechanics we 
would only have the case described by Eq. ([2]), since the (conditional) probability (density) 
to pass through one of the slits when the other is also open is just zero. 

In the literature it is common to associate the two contexts described above with either 
the duality principle (with the two slits open simultaneously, the system behaves as a wave; 
with only one, as a corpuscle) or the uncertainty principle (determining the particle position 
is the same as to determine which slit the particle passed through, while letting it pass 
through both slits open allows us to determine its momentum). However, we can also find 
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a relationship between the change of context and the problem of decoherence. As is well 

CI 

known, decoherence is nowadays the most widespread mechanism used in the literature |9J to 
explain the appearance of classical-like behaviors in quantum systems due to their dynamical 
interaction with an environment. If system plus environment are initially described by a 
separable total wave function, their coupling makes this wave function to become entangled 



10, 111) in time. When we "look" at the reduced system dynamics, this entanglement 



translates into a partial or even total loss of the system coherence. In the case of the two 
slits, this process could be seen as a sort of smooth transition from a context where both 
slits are open simultaneously, to another in which only one slit is open at a time, because 
of the gradual screening of the information which establishes that both were open when the 
particle got diffracted. In other words, as the entanglement becomes stronger, the effect is 
similar to assume that the particle track becomes more localized in time. 

The transition from one context to the other due to decoherence can be very well vi- 
sualized in terms of quantum trajectories {12I 13], since they offer a clear interpretational 
advantage by making possible to follow the system dynamics and to (intuitively) under- 
stand the underlying physics of the process at the same level that classical trajectories 
do in classical mechanics. It is well known that the standard version of quantum me- 
chanics (Copenhagen interpretation) can be reformulated in an exact fashion in terms of 
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a trajectory-based theory, namely Bohmian mechanics Jj], 

. ._ . , 

among 

many other different applications. However, and as it also happens in standard quantum 



14| . which has been widely 
used, precisely, to describe quantum interference and diffraction jisi fill 17, Q, 19] 



mechanics, Bohmian mechanics cannot be directly applied to decoherence problems with a 
arg e number of degrees of freedom because any calculation becomes numerically prohibitive 



20| | . This inconvenience can be avoided, nevertheless, by following the same philosophy 



as in standard quantum mechanics when Markovian conditions apply. In such cases, the 
effects of the environment on the system dynamics appear through some dissipative terms 
(dissipators) in the system equations of motion, as can be seen in different approaches based 
on solving either master equations for the reduced density matrix [21] or stochastic wave 



equations |22j. This allows one to devise simple, phenomeno 
the way how the coherence damping takes place 



23 



24 



25 



ogical models to understand 



26 



27 



2 a ] . Similarly, one 



can also devise simple quantum trajectory models grounded on Bohmian mechanics, such 
as the reduced quantum trajectory formalism 29J, where the environment effects appear 
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in the system equations of motion through some dissipative terms and, therefore, it is not 
necessary to deal with the dynamics of the whole system-plus-environment set of degrees of 
freedom. In this regard, we would like to mention that former analyses of decoherence by 
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3l| , where a coherent 



means of quantum trajectories were carried out by Na and Wyatt 
superposition was coupled to a harmonic bath. 

Finally, it is worth pointing out that the questions addressed here, i.e. entanglement and 
decoherence, have recently received a great deal of attention from the theoretical chemistry 



community. This arises from the feasibility of using excited vibrational states as a 



Dasis for 



quantum computing [32 



33 



35 



36 



37|, 



34j . Then problems such as entanglement dynamics 
38| . or the role of the different characteristics of the corresponding potential energy surfaces 
391 ] have been addressed in molecular states in H 2 0, S0 2 , 3 , and formaldehyde. 

In this work we explore the properties of reduced quantum trajectory formalisms in 
connection with the problems of contextuality and decoherence. In particular, we propose a 
simple trajectory model which includes also a simple mechanism of information screening. 
Although approximate, the quantum trajectories within this model are able to reproduce 
satisfactorily the projection of the "true" Bohmian trajectories onto the subspace of the 
reduced system, since it can be expected that, due to entanglement, such projections violate 
the non-crossing property of Bohmian mechanics {^(J. To make this paper self-contained, 
prior to the description of our model, we briefly analyze in Sec. [TT]the problem of decoherence 
in two-slit experiments from the standard quantum-mechanical point of view. Then, in 
Sec. IHII we describe the reduced quantum trajectory model based on the screening mechanism 
mentioned above. Numerical simulations based on this model are presented and discussed 
in Sec. IIVI A final discussion and the main conclusions derived from the present work are 
summarized in Sec. |V] 



II. DECOHERENCE IN THE DOUBLE-SLIT EXPERIMENT 

Consider a particle beam after getting diffracted by two slits. In the absence of interac- 
tions with an environment, its time-evolution can be described at any time by 

|*J 0) > = ci|Va,t> + C2[^2, t >, (3) 
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where |ci| 2 + | C2 1 2 = 1. In the coordinate representation, the element (r, r') of the associate 
density matrix, p^ = will read as 



Pi 0) (r,r') = *i 0) ( 



(4) 



and the diagonal terms giving the probability density as 

p £ (0) (r) = | Cl | 2 |-0i, £ (r)| 2 -h |c 2 | 2 |^ 2jt (r)| 2 + 2 1 Cl 1 1 c 2 1 |^ x ,, (r) | |^ 2ii (r) | cos^(r), (5) 

with fyf^ (r) = (rl^ - 1 ). In this last expression, 5 t is the space and time dependent phase-shift 
between the two partial waves. 

In order to include the effects of an environment over the system, consider now the follow- 
ing simple model (though general enough to study other interference process than two-slit 
experiments). Under the presence of such an environment, Eq. ([3]) is no longer valid to 
describe the system dynamics. Assuming that system and environment are initially decou- 
pled, the (initial) total wave function accounting for both can be expressed as a factorized 
product of the initial wave function describing each subsystem, 

|*> = |*(°)>®|^>, (6) 

with as in Eq. ([3]). If the environment acts differently on each branch of the diffracted 



beam 



2m 



271 ]. at any subsequent time we will find 

\%) = ci|^i,t) ® \E ht ) + c 2 \^ t ) <g> |£ 2 , t ), (7) 



i.e. the total wave function has become entangled. The reduced probability density associated 
with the system of interest is now obtained from (J7J) by tracing the full density matrix, 
fit = \^t){^t\i over the environment states, which leads to 

Pt= J2( E MEj,t). (8) 

3=1,2 

Note that, in the coordinate representation and for an environment constituted by N parti- 
cles, Eq. (|Sj) becomes 

Pt(r,r') = J (r, ri,r 2 , . . . r N \^ t )(^ t \r', ri, r 2 , ...r N ) dr x dr 2 ■ • -dr N . (9) 

Substituting ([7j) into Eq. (jSJ), rearranging terms and then expressing the final result in the 
reduced coordinate representation, we reach 

^(r,r') = (1 + H 2 ) l^ffeWVl^r') + 2a t cic^i,t(r)^ 2 * t (r') + c.c, (10) 

j=l,2 



where a t = (E2 t t\Ei jt ) and c.c. indicates the complex conjugate of the second term on the 
r.h.s. Finally, the (reduced) probability density resulting from Eq. ffTO]) is 

p t = (1+KI 2 ) Oci|V^r+|c 2 | V2,t| 2 + 2A t |ci||c2||^ M ||^2,t|cOS<q , (11) 

with 



M = T^T (12) 



2| 


Oit\ 




(1 + 


\oc t 


I 2 ) 



being the coherence degree 26|, \2j\ , which gives an idea of the fringe visibility of the inter- 



ference pattern. For instance, if we consider \a t \ = e */ Tc , we find 

A t = sech(£/T c ), (13) 

where r c is the coherence time, which is a function of different physical parameters (the 
system mass, temperature, etc.). In the literature one can find detailed models which a. 



low 



one to get an estimation of coherence times in different physical situations [24], |25|, |28 |. 
Nevertheless, as can be seen in Eq. fTTTl) . if r c — > 00, the decoherence process is very slow 
and the interference pattern is always well defined. On the contrary, for r c finite we observe 
an asymptotic decay of the interference pattern to a classical-like one described by 

Pt=\c l \ 2 \^,t\ 2 + \c 2 \ 2 \^ t \ 2 . (14) 

That is, we observe a smooth transition from the context where both slits were open simul- 
taneously to another one where only one slit is open at a time, not precisely because this 
was the case, but because the environment makes the information of simultaneity to become 
screened. 



III. A SIMPLE REDUCED QUANTUM TRAJECTORY MODEL FOR LOSS OF 
"WHICH- WAY" INFORMATION 

In order to visualize and describe the action of the environment over the system in terms 
of trajectories, it is convenient to express the probability current density J t in terms of the 
system density matrix, 

J t = -Im[V r pS 0) (r,r')] , (15) 

Ul r'=r 
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instead of starting from the standard Bohmian derivation in terms of \1/| ^ 40j. Thus, under 
the presence of the environment, Eq. ([15]) becomes 



h 



Jt = — Im[V r p t (r,r')] 



which satisfies the (reduced) continuity equation 

dp t 



at 



+ V3 t 



0. 



(16) 



(17) 



After Eqs. (fT5"j) and (IT7|) . one can think of the reduced probability current density as a 
transport effect of the reduced probability density through a (reduced) velocity field, v, 
according to the relation 

J t = ptv. (18) 

This field is the reduced analog of the Bohmian velocity field and we can obtain reduced 

291 ]. in analogy to Bohmian mechanics, the equation 



quantum trajectories from it by defining 
of motion 



h Im[V r p t (r, r')] 



(19) 



m Re[p 4 (r, r')] 

Due to the continuity equation (1171) and definition (1181) . the dynamics described by 
Eq. (fT9"l) leads to the correct intensity pattern when the statistics of a large particle ensemble 



is considered 



291 ]. as also happens in standard Bohmian mechanics. However, despite the 



insight on decoherence processes provided by this equation as well as its computational ad- 
vantages, a close analysis shows us that it still keeps the non-crossing property of Bohmian 
mechanics. If we consider the limit t ^> t c (a t — > 0), Eq. (fl9l) becomes 



|ci| 2 pi, f ri it + |c 2 | 2 p 2 ,tr 2 ,t 

Pcl.i 



(20) 



where ijj and pj yt are, respectively, the velocity field and probability density associated with 
the wave and 

Pc\,t = \ci\ 2 pi,t + |c 2 | 2 p2,i- (21) 
Note that, in this limit, both p c \ tt and the probability current density, 



Jci,t = Pc\,t it = \ci\ 2 pi,t h,t + \c2?p2,t h,t, 



(22) 



are properly defined. However, Eq. (|20|) still contains information on both slits, which leads 
to the non-crossing of the reduced trajectories, while the expected behavior would be just a 



crossing as will happen with the projection of the true Bohmian trajectories on the subspace 



of the reduced system 4l|. In order to observe a full transition towards a classical-like regime 
within the theoretical framework of the model here described, it is therefore necessary the 
gradual screening of such an information. 

In order to observe such a behavior, we can further proceed with our model as follows. 
If, from a trajectory point of view, a particle does not pass through one slit, we will call 
this the empty slit. Thus, in analogy to at, assume that the influence of the empty slit 
on the particle motion decreases exponentially due to the increasing entanglement with the 
environment. At the same time, the loss of information about the empty slit strengthes 
the information about the traversed one. In other words, the weight or influence of the 
empty slit decreases with time in the corresponding trajectory evolution, while that of the 
traversed slit increases. Consider that the decay of information goes exponentially. That 
is, if the particle passes through, say, slit 1, we assume that the coefficient associated with 
\ip2,t) is given by c' 2 = C2e~ t ^ Ta , where t s is the screening time, in analogy to the coherence 
time t c . Here, note that r s gives a timescale related to how fast the information provided 
by the empty slit is screened or decoupled from the particle motion due to the environment. 
On the other hand, for the coefficient of l^i,*) we choose 



c^cx^l-hPe-^VN 2 , (23) 

indicating the increasing role of the traversed slit in the quantum motion. Now, the evolution 
of the system is then described by Eq. f fl9|) . but replacing the c; by their respective time- 
dependent counterparts, c-. 

Since there are two characteristic times, t s and r c , we can define the ratio rj = t s /t c . 
Thus, if rj <C 1, the screening of the empty-slit information takes place much faster than the 
process that leads to the quenching or damping of the interference fringes. In this case, if 
the screened slit is, say, 2, Eq. ffT9l) reduces to 

r t = h,t (24) 

and the trajectories will evolve like if there was no other slit at all, i.e. like in a context 
where there is only one slit open at a time. This means that particles are allowed to cross 
the symmetry axis of the experiment because, at a given time, the momentum can have 
two different values on the same space point, as in classical mechanics. In next section we 
illustrate this processes by means of some numerical simulations. 
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FIG. 1: a) Intensity obtained from the statistics of reduced quantum trajectories (full circles) and 
standard quantum mechanics (solid line) with r c = 2.26 x 10 -2 s. b) Sample of reduced quantum 
trajectories illustrating the dynamics associated with the results shown in part (a). 

IV. NUMERICAL SIMULATIONS 



The parameters (dimensions of the two-slit assembly, particle masses and wavelengths) 
utilized in the simulations that we will present here refer to the two-slit experiment with cold 
neutrons carried out by Zeilinger et al. [3j, which we will utilize as a working model. Thus, 
in Fig. QJa), we show the excellent agreement between the statistics over reduced quantum 
trajectories (full circles) and the corresponding standard quantum mechanical calculations 
(solid line), in accordance to the continuity equation ( fTTl) and definition ( TT8l) for the reduced 
field. A sample of reduced trajectories illustrating the dynamics associated with the results 
of Fig. Ufa) is displayed in part (b) of the same figure. Here, we see that as the outgoing 
neutron beams start to interfere (at a distance of ~1 m from the two slits), some trajectories 
(mainly those closer to the symmetry axis of the experiment) start to show a conspicuous 
change of direction, which is typical in the Bohmian description of interference processes 
15 . I17I . \\\ . However, due to the interference quenching, the extent of this behavior is 



reduced in both space and time: in space because the interference effects are stronger for 
the innermost trajectories [which give rise to the central peaks in Fig. Ufa)] and in time 
because the time-of-flight of the neutrons (t/ = 2.33 x 10~ 2 s) is slightly larger than r c . 

In Fig. [2] we show the results for the extreme case of maximal decoherence, that takes 
place for r c = 0. Similarly to what we did in Fig. (TJ the statistical results obtained by 
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FIG. 2: a) Intensity obtained from reduced quantum trajectory (full circles) and standard quantum 
mechanics (solid line) calculations for a double-slit experiment with cold neutrons and r c = (no 
coherence), b) Sample of reduced quantum trajectories illustrating the dynamics associated with 
the results shown in part (a). 

means of quantum trajectories and standard quantum mechanics (a), as well as a sample 
of representative trajectories (b) are plotted. Notice that despite there is no coherence 
(in the sense that the interference terms of the reduced density matrix have been damped 
out), trajectories do not cross the symmetry axis between the two slits because they obey 
Eq. (I20p . which contains information about the two slits open simultaneously. The absence 
of interference only prevents the particles from undergoing the typical "wiggling" motion 
leading to the different diffraction peaks 

In order to illustrate how the additional mechanism leading to the screening of the in- 
formation provided by the empty slit works, we show in Fig. [3] results for several values of 
the parameter rj once the q coefficients have been replaced by the c[ ones in Eq. (|T9|) . The 
statistics of reduced-screened quantum trajectories are given in the top row and samples of 
quantum trajectories illustrating the corresponding dynamics are displayed in the bottom 
one. As can be noticed in Fig. [3ja), for small values of r\ the screening is very slow screening 
and the non-crossing dominates the system dynamics. However, as we move to higher values 
of r) (from left to right in the figure), the screening starts to increase in importance over 
the interference quenching and trajectories start to cross the axis between the two slits. In 
the extreme case, where the screening is very strong, the trajectories display a classical-like 
behavior because the information about the two slits open simultaneously initially is lost 
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FIG. 3: Top: Intensity obtained after counting of the corresponding pseudo-reduced quantum 
trajectories for: a) rj = oo, b) rj = 10, c) rj = 1, and d) rj = 0.1. The results from the counting (full 
circles) have been joined by means of B-splines (solid line) in order to facilitate their understanding. 
Bottom: Samples of trajectories illustrating the dynamics of the results shown in the upper part. 



very rapidly. 

The transition from one regime to the other one is also interesting both from the dy- 
namical and the statistical points of view. As rj increases an order of magnitude, a set of 
trajectories densely concentrates along the symmetry axis, a region where before there was 
a strong quantum force preventing the approach of trajectories. Thus, the central intensity 
maximum is remarkably enhanced, as can be seen in Fig. [3]^b), when compared with the 
same peak in Fig. [3]^a). This enhancement occurs at the expense of the intensity of the other 
maxima. That is, there is a transfer of trajectories from the outermost to the innermost 
diffraction channels. Dynamically, once the strong boundary imposed by the quantum po- 
tential along the symmetry axis (which is stronger than along the directions separating the 
different diffraction channels [3]) is suppressed, the "quantum pressure" 18, 19] pushes the 



trajectories towards the region covered by the empty slit, and favors their transfer from one 
diffraction channel to the nearby one. 

As r) is further increased, quantum trajectories are also able to penetrate more across the 
symmetry axis, as it is apparent in Fig. [3]^c). The availability of a wider accessible region 
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in the other side of the symmetry axis makes the concentration of trajectories along this 
direction to decrease, and they distribute more homogeneously. This causes a remarkable 
decrease of the coherence degree in the measured intensity, although the central peaks are 
still relatively intense. This trend continues until the interference pattern completely dis- 
appears when the decoupling is maximum [see Fig. [3(d)]. In this case, the trajectories are 
unaffected by the presence of the other slit, i.e. they display a totally classical-like behav- 
ior. Moreover, a full transition from a dynamics characterized by single-valuedness of the 
momentum to another where it is bi- valued is observed, taking place this process within the 
system subspace. 



V. FINAL DISCUSSION AND CONCLUSIONS 



Trajectory-based approaches have received much attention in the last 



years as a potential 



tool to handle and study high- dimensional complex quantum systems [421 ] . Nowadays the 
design of new numerical tools highly relies on this kind of formalisms rather than using 
other approaches based on the time-dependent Schrodinger equation [20|. In the particu- 
lar case of Bohmian mechanics, this computational power combines with its capability to 
provide causal interpretations to quantum phenomena within a purely quantum-mechanical 
framework, unlike other semiclassical or classical approaches. These two interesting features 
have brought Bohmian mechanics from the field of the Foundations of Physics to the most 
applied fields in physics and chemical physics, attracting the interest of many different com- 
munities within. Similarly, the model proposed here to study decoherence and the process 
of going from one context to another one in quantum mechanics provides an understandable 
and intuitive insight into dynamics involved. In particular, regarding the mechanism that 
leads a linear theory with single-valued momenta to an "apparent" non-linear one, where 
the momentum can be multi-valued. 

Here we have dealt with the problem of the damping or quenching of the interference 
fringes produced by decoherence in a two-slit experiment under the presence of an environ- 
ment, which yields as a result a classical-like pattern. This effect is, to some extent, similar 
to consider the transition from a context where both slits are open simultaneously to another 
one where only one slit is open at a time and, therefore, the particle track is localized (i.e. 
the particle has passed through one slit or the other). In order to elucidate and understand 
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decoherence without taking into account explicitly the dynamics of the environment degrees 
of freedom, we have considered some simple reduced quantum-trajectory models. Though 
limited because of their simplicity, one can infer from them that, unless there is a sort of 
screening of the empty-slit information, the patterns affected by decoherence can be well 
reproduced, but keeping still a sort of internal coherence which makes the trajectories to 
avoid their crossing, as in standard Bohmian mechanics. Only when the empty-slit informa- 
tion is gradually screened, the trajectories start to cross, as one would expect when using 
Bohmian mechanics and projecting the true 3(A^+ l)-dimensional quantum trajectories onto 
the subspace of the reduced system. Thus, though simple and approximate — more refined 
and precise models, which are out of the scope of this work, but that are being currently de- 
veloped in order to go beyond the analysis presented here — , this model allows us to study 
the properties of decoherence at a relatively cheap computational cost and providing, at 
the same time, a physical insight on the way how the contextual dependence and nonlocal 
correlations are gradually suppressed. 
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